ITEP/TH-68/09 

New and Old Results in Resultant Theory 

A.Morozov * and Sh.Shakirov ''' 
ABSTRACT 

Resultants are getting increasingly important in modern theoretical physics: they appear whenever one deals with non-linear 
(polynomial) equations, with non-quadratic forms or with non-Gaussian integrals. Being a subject of more than three-hundred-year 
research, resultants are of course rather well studied: a lot of explicit formulas, beautiful properties and intriguing relationships are 
' known in this field. We present a brief overview of these results, including both recent and already classical. Emphasis is made on 
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explicit formulas for resultants, which could be practically useful in a future physics research. 
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1 Introduction 



This paper is devoted to resultants natural and far-going generalization of determinants from conventional 
linear to generic non-linear systems of algebraic equations 

fi{xi,...,Xn) = 

f2{xi,...,Xn) = 

< 

where fi, - ■ ■ ,fn are homogeneous^ polynomials of degrees n, . . . , r„. This system of n polynomial equations 
in n variables is over-defined: in general position it does not have non-zero solutions {xi, . . . ,Xn) at all. A 
non- vanishing solution exists if and only if a single algebraic constraint is satisfied, 

R{fl,f2,...,fn}=0 (1) 

where R is an irreducible polynomial function called resultant, depending on coefficients of the non-linear system 
under consideration. Clearly, for linear equations resultant is nothing but the determinant - a familiar object 
of linear algebra. For non- linear equations, resultant provides a natural generalization of the determinant and 
therefore can be considered as a central object of non-linear algebra, [1]. 

Historically, the foundations of resultant theory were laid in the 19th century by Cayley, Sylvester and 
Bezout [2]. In this period mostly the case of n = 2 variables was considered, with some rare exceptions. The 
theory of resultants for n > 2 and closely related objects, namely discriminants and hyperdeterminants, was 
significantly developed by Gelfand, Kapranov and Zelevinsky in their seminal book [3]. 

Let us emphasize, that determinant has a lot of applications in physics. Resultant is expected to have at 
least as many, though only a few are established yet. This is partly because today's physics is dominated by 
methods and ideas of linear algebra. These methods are not always successful in dealing with problems, which 
require an essentially non-linear treatment. We expect that resultants and related algebraic objects will cure 
this problem and become an important part of the future non-linear physics and for entire string theory program 
in the sence of [4]. To illustrate this idea, let us briefly describe the relevance of resultants to non-Gaussian 
integrals - a topic of undisputable value in physics. 

It is commonly accepted, that quantum physics can be formulated in the language of functional integrals. 
It is hard to question the importance of this approach: since the early works of Wigner, Dirac and Feynman, 
functional integrals have been a source of new insights and puzzles. Originated in quantum mechanics, they 
have spread fastly in quantum field theory and finally developed into a widely used tool of modern theoretical 
physics. Prom the point of view of functional integration, any quantum theory deals with integrals of a form 

=Jdx fix) e-^(-) (2) 

^In this paper we assume homogeneity of the polynomial equations to simplify our considerations; actually, the homogeneity 
condition can be relaxed and the whole theory can be formulated for non-homogeneous polynomials as well. 
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called correlators. The integration variable x belongs to a linear space (which is often infinite-dimensional), and 
the choice of particular action S{x) specifies a particular theory or a model. In the simplest models, S{x) is 
merely a quadratic form, which in some particular basis can be written as 



S{x) — ^ ^ SijXiXj (3) 



If this is the case, the integral (2) is called Gaussian. In a sence, Gaussian integrals constitute a foundation 
of quantum theory: they describe free (non-interacting) particles and fields. These systems are relatively 
simple, and so are the Gaussian integrals. For example, the simplest two-dimensional Gaussian integral can be 
evaluated, say, by diagonalisation, and equals 

J J \/}r — Aac 

where the expression 6^ — 4ac can be recognized as a determinant of the quadratic form ax'^ + hxy + cy^ . Other 
Gaussian integrals, with more variables, with non-homogeneous action or with non-trivial insertion /(</>), are 
equally easy and similarly related to determinants - central objects of linear algebra, which describe consistency 
of linear systems of equations. Remarkably, the study of low-dimensional non-Gaussian integrals of degree three 

SijkXiXjXk 

(5) 

and higher, reveals their relation to resultants: for example, a two-dimensional non-Gaussian integral equals 



hi 



where the expression TJa^iP — b^c^ — ISabcd + Aac? + Ab^d can be recognized as a discriminant of the cubic 
form ax^ + bx^y + cxy^ -\- dy^, given by resultant of its derivatives: 

ii| Sax'^ + 2bxy + cy"^, bx"^ + 2cxy + My^ | = 270^^^ - b'^c^ - ISabcd + 4ac^ + Ab^d (7) 



In this way, resultants naturally appear in the study of low-dimensional non-Gaussian integrals. This connection 
between functional integration and non-linear algebra can be traced much further: section 5 of present paper 
provides a brief summary of this interesting direction of research. See also [5] for more details. 

In this paper we present a selection of topics in resultant theory in a more-or-less explicit manner, which is 
more convenient for a reader who needs to perform a practical calculation. This also makes our paper accessible 
to non-experts in the field. Most of statements are explained rather than proved, for some proofs see [3] . 
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2 Basic Objects 
Homogeneous polynomial 

Basic objects of study in non-linear algebra are generic finite-dimensional tensors [1]. However, especially 
important in applications are symmetric tensors or, what is the same, homogeneous polynomials of degree r in 
n variables. In present paper, we restrict consideration to this case. It should be emphasized, that restriction 
to symmetric tensors is by no means an oversimplification: most of the essential features of non-linear algebra 
are present at this level. There arc two different notations for homogeneous polynomials, which can be useful 
under different circumstances: tensor notation 

n 

ii,i2,---,v=l 

and monomial notation 

To distinguish between these notations, we denote coefficients by capital and small letters, respectively. A 
homogeneous polynomial of degree r in n variables has C^_|_^_i independent coefficients. 

Resultant 

A system of n polynomial equations, homogeneous of degrees n , . . . , r„ in variables xi,. . . ,Xn 

/i(xi, . . . ,a;„) = 

f2{Xl,- ..,Xn)=0 

< 

fn (^1 5 • • • 5 ) ~ 

has non-zero solutions, if and only if its coefficients satisfy one algebraic constraint: 

Rri,...,r^{fl, ■ ■ ■ , fn} = 

The left hand side of this algebraic constraint is called a resultant of this system of equations. Resultant always 
exists and is unique up to rescaling. Resultant is a homogeneous polynomial of degree 

deg-Rr-i,...,r„ =ri...r„ ( — + ... + — ] 

in the coefficients of the system. Moreover, it is homogeneous in coefficients of each equation /j separately, and 
has degree ri . . . r^-iri+i . . . r„ in these coefficients. When all equations are linear, i.e, when = 1, resultant 
is equal to the determinant of this system: -Rii...i{/i, /2, • • • , fn} = /i A /2 A . . . A /„. Resultant is defined 
only when the number of homogeneous equations matches the number of variables: when these numbers differ, 
various more complicated quantities should be used to describe the existence and/or degeneracy of solutions. 
Some of these quantities are known as higher resultants or subresultants [6], some as complanarts [7]. 
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Discriminant 

A single homogeneous polynomial S{xi, . . . ,Xn) of degree r in n variables gives rise to a system of equations 

' dS{xi,...,Xn) 



dxi 



= 



dS{xi,. ..,x„) 

dXn 



which are all homogeneous of degree r— 1 in variables xi,. . . ,Xn- Resultant of this system is called a discriminant 
of S and is denoted as Z)„|r(iS'). Discriminant is a homogeneous polynomial of degree 



deg2)„|,(5) = n(r-l) 



n-l 



in the coefficients of S, what follows from the analogous formula for the generic resultant. When 5 is a quadratic 
form, i.e, when r = 2, discriminant is equal to its determinant: £>„|2(<S') = detS'. 

Invariant 

A Lie group of linear transformations 

/ Xi \ ( Gil Gi2 ... Gin \ / xi \ 



\ X„ ) \ G„i Gi2 ... Gnn J \ J 

with detG ^ is called GL{n). Its subgroup with detG = 1 is called SL{n). In tensor notations, they act 

Thus they act on coefficients of homogeneous polynomials S{xi, . . . , Xn) of degree r in n variables by the rule 

Si-i^...^i^ <—> Gi-ij-i . . . Gi^j^Sj-^^,,,j^ 

Polynomials I{S) which map to theirselves under SL{n) transformations are called SL{n) invariants of S. In 
particular, discriminant Dn\r{S) is a SL{n) invariant of S. However, there are many more invariants which 
have no obvious relation to the discriminant. The number of functionally independent invariants of S is equal 
to the number of independent coefficients of S minus the dimension of SL{n), that is 

# functionally independent invariants of S = C^_|_^_i — + 1 



Therefore, invariants of given degree form a finite-dimensional linear space. All invariants form a ring. 



5 



3 Calculation of resultant 



3.1 Formulas of Sylvester type 



To find the resultant of n polynomial equations, homogeneous of degree r in variables xi,. . . ,x. 



"n 



Xn) = 



fliXi 



Xn) = 



< 



^ fn{x\ 5 • • • 5 Xn) — 

it is possible to consider another system, multiplied by various monomials: 

consider x^ . . . x^ fi =0 instead of just fi = 

In the monomial basis, coefiicients of this system can be viewed as a rectangular nC^_,_„_2 x C^^^'_|_„_2 matrix, 
where g = si + . . . + s„ is the multiplier's degree. This matrix is called the (generalized) Sylvester matrix. 

3.1.1 Resultant as a determinant of Sylvester matrix 

In some cases it is possible to choose q in such a way, that Sylvester matrix is square, i.c, nC^_^_jj^_i = Cl^Xq+n-i- 
In these cases, resultant is equal to the determinant of Sylvester matrix. Looking closer at this numeric relation, 
it is easy to see that the only cases when it is possible are n = any, r = 1 and n = 2, r = any. In the first 
case, it suffices to choose q = and Sylvester matrix equals the matrix of the system. This is the case of linear 
algebra. The second case is less trivial: it is necessary to choose q = r—1 to make the Sylvester matrix square. 
For example, in the case n = 2, r = 2 the system of equations has a form 

{fnxf + /i2a;ia;2 + /22a;2 = 
giixi + 5i2a:;ia;2 + 922x1 = 
Multiplying the equations of this system by all monomials of degree q = 1, that is, by Xi and X2, we obtain 

flixl + fl2xlx2 + f22Xlxl + 0x1 = 
0x1 + fnXiX2 + f\2X\x\ + h2x\ = 



9nXi + 9i2x\x2 + 922X1x1 + Oa;^ = 

Oxf + giixlx2 + 9i2Xixl + 922x1 = 

These four polynomials are linear combinations of four monomials, i.e, give rise to a square Sylvester matrix: 



< 



^2|2{/,5} = det 



/ hi fl2 f22 \ 
/ll /l2 /22 



(8) 



ffll 512 922 







\ gii gi2 922 I 
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Formulas, which express resultant as a determinant of some auxiliary matrix, are usually called determinantal. 
In the next case n = 2, r = 3 the system of equations has a form 



fllixl + fll2xlx2 + hllX\x\ + !l1ix\ = 

g\\\x\ + gwixix-i + gv2,ix\^x\ + 52222:2 = 
Multiplying the equations of this system by all monomials of degree q = 2, that is, by xf, X1X2 and X2, we get 

/111a;? + hl2x\x2 + fl22xlxl + f222xlxl + Oxixf + Ox^ = 

Oxf + /iiia;fa;2 + fii2xlxl + fi22xlxl + /222a;ia;| + Oa;| = 
Qx\ + Qx\x2 + hiix\xl + hi2xlxl + /i22a;ixt + /222a;i = 

51113^1 + g\\2x\x2 + g\22x\xl + g222xlxl + Oxixf + Ox^ = 

^xl + giiixf a;2 + gii2x\xl + gi22x\xl + 52222:12;^ + 0x1=0 
Ox\ + Ox\x2 + 51112^1 2^2 + gii2xlxl + gi22Xixl + 5222a;2 = 

These four polynomials are linear combinations of four monomials, i.e, give rise to a square Sylvester matrix: 



R2\i{f,g} = det 



/ /ill 


/112 


/122 


/222 





\ 





/ill 


/112 


,/l22 


./222 











/ill 


,/ll2 


./122 


/222 


fflll 


5112 


.9122 


.9222 











.9111 


.9112 


.9122 


.9222 





V 





5111 


5112 


5122 


5222 / 



(9) 



Note, that the matrix has a typical ladder form. The beauty and simplicity made this formula and corresponding 
resultant i?2|r widely known. In our days, Sylvester method is the main apparatus used to calculate two- 
dimensional resultants and it is the only piece of resultant theory included in undergraduate text books and 
computer programs such as Mathcmatica and Maple. However, its generalization to n > 2 is not straightforward, 
since for n > 2 and r > 2 Sylvester matrix is not square. 

It turns out that, despite Sylvester matrix in higher dimensions is rectangular, it is still related closely to 
the resultant. Namely, top-dimensional minors of the Sylvester matrix are always divisible by the resultant, 
if their size is bigger, than resultant's degree (i.e, for big enough q). Various generalizations of the Sylvester 
method explicitly describe the other, non-rcsultantal, factor. Therefore, such generalized Sylvester formulas 
always express the resultant as a ratio of two polynomial quantities. This is certainly a drawback, at least from 
the computational point of view. 

Sylvester's method posesses at least two different generalizations for n > 2: one homological, which leads 
to the theory of so-called Koszul complexes, another more combinatorial, which leads to the construction of 
so-called Macaulay matrices. 
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3.1.2 Resultant as a determinant of Koszul complex 

In construction of the Koszul complex, one complements the commuting variables 



with anti-commuting variables 

^1, 62, . . . ,6n ' 616 j + 6j6i = 

and considers polynomials, depending both on xi, . . . ,Xn and 6i, . . . ,9n- Denote through Cl{p, q) the space of 
such polynomials of degree p in a;- variables and degree q in ^-variables. The degree q can not be greater than 
n, since 9 are anti-commuting variables. Dimensions of these spaces are 

dimO(p,g) = C^+„_iC« (10) 

Koszul differential, built from /, is a linear operator which acts on these spaces by the rule 

d d 

d = fl{xi, ^n)-g^ + ... + fn{xi, . . ■^^n)-Q^ (11) 

and is automatically nilpotent: 

d d f d d d d \ 

It sends 

d: n{p,q)^n(p + r,q-l) 

giving rise to the following Koszul complex: 

VL{p,q) ^ n{p + r,q - I) ^ nip + 2r,q - 2) 4 ...4 r!(0,i?) 

Thus, for one and the same operator d there are many different Koszul complexes, depending on a single integer 
parameter R - the a;-degree of the rightmost space. Using the formula (10) for dimensions, it is easy to write 
down all of them. For example, the tower of Koszul complexes for the case 3|2 is 
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R 



Spaces 



Dimensions 



Euler characteristic 



312 



n{o,o) 

1 ri(i,o) 

2 0(0,1) ^ Q(2,0) 

3 0(1,1)^0(3,0) 

4 0(0,2) -> 0(2, 1) ^ 0(4,0) 

5 0(1,2)^0(3,1)^0(5,0) 

6 0(0,3)^0(2,2)^0(4,1)^0(6,0) 

7 0(1,3)^0(3,2)^0(5,1)^0(7,0) 



1 
3 

3^6 

9^ 10 
^18^ 
^30^ 

18^ 45 
30^ 63 



15 
21 



28 
36 



where Euler characteristic is the alternating combination of dimensions. The rightmost differential 



d: 0(i?-r,l) A O(ii,0) 



acts on the corresponding linear spaces as 



(^XiX2 ■ ■ ■ x]^^ • Oj I— > (^XiX2 ■ ■ ■ a^n"^ • fj {xi, . . . , Xn) 



and is represented by nC^Zr+n-i ^ ^R+n-i matrix. It is easy to see, that this is exactly the rectangular 
Sylvester matrix. Therefore, Koszul complex in a sence complements the Sylvester matrix with many other 
matrices, which carry additional information. Resultant is equal to certain characteristic of this complex 



-Rn|r{/li • • • ) fn} = DET (^1,(^2, . . . , tip-l) 



(12) 



which is also called determinant, not to be confused with ordinary determinant of a linear map. Here di are 
the rectangular matrices, representing particular terms of the complex. The " determinantal" terminology was 
introduced by Cayley, who first studied DET (di,d2, . . . ,rfp_i), and is rather imusual for modern literature, 
where DET {di,d2, ■ ■ ■ , dp-i) is called a (Reidemeister) torsion of the complex di, . . . , dp-i. Since the quantity 
itself is a rather natural generalisation of determinants from linear maps to complexes, we prefer to use the 
historical name " determinant of a complex" . 

Calculation of determinants of complexes is described in some detail in [3, 8]. Given a complex 



d2 



(13) 



with vanishing Euler characteristic 



X = dim L\ — dim L2 + dim L3 — . . . + (—1)^ dim Lp = Q 



one needs to select a basis in each linear space Lj and label the basis vectors by elements of a set 



9 



{1,2,. ..,k} 



After the basis is chosen, the complex is represented by a sequence of rectangular matrices. The U x Zj+i matrix 
di has k rows and k+i coloums. To calculate the determinant of this complex, one needs to select arbitrary 
subsets cTj C {l, 2, . . . , Zj} which consist of 



|a,|=^(-ir+'=+iZ, 



fe=i 



elements, i.e. 



CTl = 0, <T2 = ^1, CT3 = Z2 — ^1, kr4 = Z3 — ^2 + ^1 



and define conjugate subsets CTj = {l, 2, . . . , ij}/(7i which consist of 



elements. Finally, one needs to calculate minors 



Mi = the minor of di, which occupies rows crj and coloumns c7,+i 



The determinant of complex is then given by 



p-i 



DET {di,d2,...,dp-i) = l[{Mi) 



{-!)" 



Mp_2Mp_4 . 



(14) 



In fact, this answer docs not depend on the choice of o",; [8] and thus gives a practical way to calculate resultants 
of polynomial systems. Given a system of equations, one straightforwardly constructs the diifcrcntial (11) and 
the linear spaces (10), calculates the matrices and their minors according to the above procedure and, finally, 
uses the formula (14) to find the resultant. 

The answer, obtained in this way, is excessively complicated in the following two senses: (i) the ratio of 
various determinants at the r.h.s. is actually reduced to a polynomial, but this cancelation is not explicit in 
(12); (ii) the r.h.s. is actually independent of the choice of Koszul complex, but this independence is not explicit 
in (12). Still, Koszul complex provides an explicit tool of calculation of resultants, which is reasonably useful 
for particular low n and r. 
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To illustrate the Koszul complex method, let us consider in detail the case n = 3, r = 2, a system of three 
quadratic equations in three unknowns. Let us denote their coefficients by a, b, c instead of usual /, g, h: 

fi{xi,X2,X3) = aiixl + ai2Xia;2 + 013X1X3 + a22x| + 023X2X3 + a33x| 

/2(xi, X2, X3) = biixl + &12X1X2 + &i3a;ia;3 + ^223:2 + b23X2X3 + hsxl 

f3{xi,X2,X3) = CiiX? + C12X1X2 + C13X1X3 + C22a:;i + C23a;2a;3 + C33X^ 



By definition, Koszul differential is given by 



Let us describe in detail the simplest case R — 4 and next-to-thc-simplest case R — 5. The complex with R = A, 
is the first determinant-possessing complex (i.e. has a vanishing Eulcr characteristics): 



n{0,2) n{2,l) ^ 17(4,0) 



with dimensions 3-^18— >15. If we select a basis in il(0, 2) as 



a basis in f2(2, 1) as 



{0203, 0103, 0102} 



{Xi03,X^02, Xi6'i, XiX26'3, XiX26'2, XiX26'i, XiX36'3, XiX36'2, 



XiX36'i, X26'3, X26'2, .X26'i, X2X36'3, X2X36'2, X2X36'i , X36'3 , X36'2 , a;36'l} 



and, finally, a basis in f7(4,0) as 



r43 3 2 2 2 22 3 2 2 343 2 2 34-1 

^Xi, X1X2, X^^X^^ X-yX'2-j X1X2X3, '^l'^3; X\X'2^ X1X2X3, X1X2X3, X1X3, X2 , -^2^3' ^2^3, '^2X31 X'^j 



then Koszul differential is represented by a pair of matrices: one 3 x 18 



■^11 ''11 ^ ^12 ''12 ^ ^13 '"13 ^ ^2 2 '"2 2 ^ '^2 3 ^^2 3 ^ '^3 3 '"3 3 



—b-.-. — T — 6. r, — o — o —a.-, .-, — i.-, .-, 



■"2, 2 "^2, 2 ^ "2, 3 "2, 3 



2,3 ^"3, 3 



3. 3 
3, 3 "-^3, 3. 
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and another 18 x 15 



=1.1 


= 1, 2 


= 1,3 


=2.2 


=2, 3 


=3.3 





























-*1.1 


-*1,2 


-*1.3 


-*2.2 


-*2,3 


-*3,3 





























"1.1 


■=1,2 


"1,3 


"2.2 


"2,3 


"3.3 
































= 1, 1 





= 1.2 


= 1,3 





=2, 2 


=2, 3 


=3.3 























-*1,1 





-*1.2 


-i, 3 





-*2,2 


-*2,3 


-*3.3 























"1, 1 





"1.2 


"1,3 





"2, 2 


"2,3 


"3.3 


























=1,1 





= 1, 2 


= 1.3 





=2, 2 


=2.3 


=3, 3 























-*1,1 





-*1,2 


-*1,3 





-*2,2 


-*2.3 


-*3,3 























"1,1 





"1,2 


"1.3 





"2, 2 


"2.3 


"3,3 


























= 1.1 








= 1,2 


= 1, 3 








=2,2 


=2, 3 


=3. 3 

















-*1.1 








-*1,2 


-*1,3 










-*2,3 


-*3,3 

















"1.1 








"1,2 


"1,3 








"2.2 


"2,3 


"3,3 




















=1,1 








= 1, 2 


= 1.3 








=2, 2 


=2, 3 


=3.3 

















-*1,1 








-*1,2 


-*1.3 








-*2,2 


-*2,3 


-*3.3 

















"1,1 








"1,2 


"1.3 








"2, 2 


"2,3 


"3.3 




















=1.1 








= 1.2 


= 1,3 








=2. 2 


=2.3 


=3,3 

















-*1,1 








-*1.2 


-i, 3 








-*2,2 


-*2.3 


-*3,3 

















"1.1 








"1.2 


"1,3 








"2,2 


"2.3 


"3,3. 



By selecting some 3 columns in di and complementary 15 rows in c?2, we obtain the desired resultant 



"1. 1 

"1.2 
"1,3 
"2,2 
"2,3 
"3,3 



'1,2 
= 1,3 




"1.1 

"1,2 " 

"1,3 "^1,2 

-^1,3 






'2, 2 


"2,2 














'2, 3 


"2.3 


=2, 2 


-^2.2 


"2.2 





=3.3 


"3.3 


=2, 3 


-*2.3 


"2.3 











=3,3 


-*3,3 


"3.3 





































































































































































1, 1 


-'■1,1 


"1,1 


























=1, 1 


-*1,1 


"1,1 


























=1, 1 




1, 2 


-^1,2 


"1,2 

















1, 3 




"1.3 


= 1,2 


-^1.2 


"1.2 

















= 1,3 


-*1,3 


"1.3 


= 1, 2 


-*1, 




















= 1,3 


-*1, 



=2, 2 
^2. 3 
=3,3 




'■2.2 
*2,3 
■'■3,3 




'2.2 
'2,3 
'3,3 



2, 2 


-^■2,2 


"2,2 











2, 3 


-'■2,3 


"2,3 


'2, 2 


-^■2,2 


"2,2 


3, 3 


-'3,3 


"3,3 


'2, 3 


-'2,3 


"2,3 











'3, 3 


-*3,3 


"3,3. 



*1, 1 "^1,1 "^1,2 
-<j] J 
-a, , -a, T 



This particular formula corresponds to the choice of columns 1, 2 and 5 in the first matrix, but the resultant, 
of course, does not depend on this choice: any other three columns will do. One only needs to care that 
determinant in the denominator does not vanish. If it is non-zero, then the upper minor is divisible by the lower 
minor, and the fraction is equal (up to sign) to the resultant. 
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Alternative complex is the second complex with vanishing x, that corresponds to i? = 5: 

n{i,2) ^ rj(3,i) ^ 17(5,0) 

It has dimensions 9 ^ 30 ^ 21. If we select a basis in fi(l, 2) as 

{Xi9203, Xi9i0z, X16162, X292O3, X2O16Z, X29l92,Xz6203,Xz0i0z,Xz6i62} 

a basis in il(3, 1) as 

{x\Oz, xl92,xlei,xlx203, xlx292,x\x20i,x\x:iB3, xlx392, xfx^Bi, 

Xixl93,xixl92, xixl9i, X1X2X393, X1X2X392, X1X2X391, xixl93, xixl92, xixl9i, x^Bj, xl92, xl9i, 
x\x393, xlx392,xlx39i,X2xl93,x2xl92, X2xl9i,xl93, xl92,xl9i 

XiX39i,xl93,xl92, xl9i,X2X393,X2X392,X2X39i,xl93,xl92,xl9i} 

and, finally, a basis in il(5, 0) as 

/54 4 323 322322 2 223 

\Xi, XiX2t X1X3, X1X2, X1X2X3, X^X^, X1X2, X1X2X3, XiX2X^, X^X^, 

X1X2, X1X2X3, XiX2X'j, XiX2X^, XiX^, X2, X2X3, X2X^, X2X^^ X2X^^ Xg} 

then Koszul differential is represented by a pair of matrices: one 9 x 30 

^1.1 " ''1.1 ^1.3 " ■'l.i ^1.5 " ^1,1 ^.3 ° ■'3,3 ° ^3.5 ^.5 ^ '5.5 " 

^1 J ^1^1 ^1^2 ^1^3 ^1^3 Q 

Q Si^i Q Q " ^3,3 

^1- 1 Q ^1^1 fj^i ^1^2 jj^j Q 

^"■'ll^ll^ " ^ 0^;j2^13^^1i^l5^ ^ ^ " 

^^ll'll^ ^ ^ OSj^trj^OS^jtrj^O 

O^jjjO i^j O^i^jO 'iS^^'li^ '15 " 








'2,1 


Q 


^.3 


"3,3 


(1 








'1,3 







J 
















^ J 







%2 


Q 




■"2,3 


















"2,^ " ■^3,3 








-^,2 


^,2 







^,3 " ^.3 ^.3. 



and another 30 x 21 



-^,2 -^2,1 



°2,2 ^,3 ^,3 



-*2,1 ->^,i ° 



° -^,2 

" ^,3 ^,3 ^.3 



° ^,2 ^.i 
" °S,S °S,3 °3,3 



f ^,3 ^.3 ^.3 " 
II fl2^3 flj^s II 



^.2 ^.5 
" ■^^S ■=3,5 " 



f ^.2 ^.3 ^.3 
II -^,1 -^.,1 



" ^.i ^.3 ^.3 

" ■=s,s ■=s,s 
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By selecting some 9 columns in di and complementary 21 rows in d2, we obtain the desired resultant 



"-^1,1 
































































^1,2 


'1,1 


=1,1 


































































^1,1 


=1,1 




















































^,2 


'1,2 


"1,2 








"^1,1 ^1,1 
















































"^1.3 


"1,3 


^1,2 


=1,2 








"1.1 


=1,1 
















































-^1,3 


=1,3 














^1,1 


=1,1 





































'2,2 


=3,2 








"1,2 ^1,2 














^1,1 


































'2,1 


"2,3 


^,2 


=2,2 




3 ^1,3 


"1,2 


=1,2 











"1,1 


^1,1 


=1,1 

























'3,1 


"3,3 




=2,3 








"1,3 


=1,3 


^1,2 


=1,2 
















1, 1 


^1,1 


=1,1 




















-^,3 


=3,3 














^1,3 


=1,3 

























"1,1 


^1,1 


=1,1 

















'2, 


2 ^,2 














^1,2 














































'2, 


3 ^,3 


"2,2 


=2,2 








^1,3 


"1,2 


^1,2 


=1,2 







































3 ^,3 


"2,3 


=2,3 


^,2 


=2,2 





"1,3 


^1,3 


=1,3 


t: 




^1,2 


=1,2 
































"3,3 


=3,3 


^,3 


=2,3 














t: 


1,3 


^1,3 


=1,3 


"1,2 


^1,2 


=1,2 





























^,3 


=3,3 

























"1,3 


^1,3 


=1,3 



































2 
































































3 




2 


"2 2 























































-h,3 


"2 5 


2 


=2 3 


"2,2 


2 


=2 2 















































"3 3 


-h,3 


=3 3 


"2,3 


3 


"2 3 




2 


"2 2 















































"3,3 


^,3 


=3 3 




3 


"2 3 


























































"3 3 


-h,3 


=3 3 














■*1,1 





"1,2 


*1,3 





"2,3 


*3,3 


































^1,1 


1, 1 





^1,3 


"2,2 





^,3 





































*1,1 


^1,2 





^.2 


^,3 











































"1,1 








"1,3 





*2,2 











































"1,2 








^,2 


"2,; 


































^1,1 





^1,2 


^1,3 








^, 


2 



































*1,1 





"1,2 


*1,3 











































^1,1 








^1,3 

















































-=1,2 


























This particular formula corresponds to the choice of columns 1, 3, 5, 7, 12, 14, 16, 19 and 21 in the first matrix. 
Of course, any other choice will give the the same (up to a sign) result, if only the minor in denominator is 
non-zero. Direct calculation (division of one minor by another) should demonstrate, that the two complexes 
3 ^ 18 — > 15 and 9 ^ 30 ^ 21 give one and the same expression for the resultant. It should also be reproduced 
by determinants of all other complexes from the 3|2 list with i? > 3. 

3.1.3 Macaulay's formula 

The second, rather combinatorial, generalization of the Sylvester method is due to Macaulay. This method allows 
to express the resultant as a ratio of just two determinants (for Koszul complexes (14) both the numerator and 
denominator are products of several determinants). For this reason, Macaulay formula is computationally more 
economical, than eq. (14). We do not go into details here, they can be found in [10]. 
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3.2 Formulas of Bezout type 

Consider a system of n polynomial equations, homogeneous of degree r in variables xi, . . . ,x„: 

fi{xi,...,Xn) = 

f2{xi,...,Xn) = 
^ fnip^l •)•••■) ^n) ~ 

where 



fiip^l ) • • • ; ^n) 



^ ^ fi\aia2---an'^l "^2 • • • 



ai+a2 + ---+a„=r 



In the monomial basis, coefBcients of this system can be viewed as forming an n x C^_^_^_i rectangular matrix. 
The set of all top-dimensional minors of this matrix is naturally a tensor with n indices: 



= nx n minor of /, situated in coloumns 



Resultant R{fi, ■ ■ ■ , /«} is expected to be a homogeneous polynomial of variables Mi^,,,i^ of degree r"~"^. Such 
expressions are called formulas of Bezout type. The simplest example of Bezout formulas is the case n = 2, r = 1, 
where the system of equations is just a simple linear system: 

flXl + f2X2 = 

gixi + g2X2 = 



and the only top-dimensional minor is 



Mi2 = hg2 - f2gi 



In this case, resultant is just equal to this minor: 



R2\i{f,g} = Mi2 



In the next-to-simplest case n = 2, r = 2 the system of equations has a form 

fnXi + fl2XlX2 + f22X2 = 

gnxj + gi2XiX2 + g22x'i = 
and there are three top-dimensional minors: 



15 



Afl2 = /llfi'12 - /1251I 
Mi3 = /11522 — f229U 
M23 = f 12922 — /22912 

In this case, resultant is a homogeneous polynomial of degree 2 in minors , which is moreover a determinant 

M12 Mis 



R2\2{f,9} = M12M23 - Mi^3 = det 



Mi3 M23 



For n = 2, r = 3 the system of equations has a form 

+ f\\2x\xi + /i22a;ix| + /222a;2 = 

g\\\x\ + g\\2x\x'2, + 9V2,2X\X\ + 5f222a;2 = 

and there are six top-dimensional minors: 

Afl2 = /lllS'112 — /ll2S'lll -^13 = /lllS'122 — /l225'lll 
Afl4 = /lllfl'222 — /222fl'lll -^23 = /1125122 — /l22fl'112 
-^24 = /1125222 — /2225112 -M34 = /1225222 — /2225122 

In this case, resultant is a homogeneous polynomial of degree 3 in minors, which is again a determinant 

^2|3{/, 5} = M12M34M14 + M12M34M23 - M12MI4 - M23M34 + 2M13M14M24 - Mf4 - Ml^M^-i (15) 



-R2|3{/,6'} = det 



/ M12 Mi3 Mi4 \ 

Mi3 M23 + Mi4 M23 

Mi4 M24 M34 



V 



(16) 



Note, that six variables Mjj are not independent: they are subject to one constraint 



M12M34 - M13M24 + M14M23 = 



(17) 



which is well-known as Plucker relation. For this reason, the polynomial expression is actually not unique: there 
are many different polynomial expressions for resultant through minors, whose difference is proportional to the 
Plucker relation. Generally (for higher n and r) there are several Plucker relations, not a single one. Variables 
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Mi^,,,i^ are sometimes called Plucker coordinates. Bezout formulas, which express resultant through Plucker 
coordinates, exist in higher dimensions: for example, for n = 3, r = 3 we have 

f{xijX2,X3) = fiixl + /i2a;ia;2 + fi3.xiXs, + f22xl + f23.X2Xz + fzs.xl = 
< g{xi,X2,xz) = giixl +5i2a;ia;2 +gi3a;ia;3 + g22xl +5232:22:3 +5332:1 = 
h{xi,X2,X3) = hiixf + /ii2a;iX2 + hisxixs + h22xl + h23X2X2, + h^d,x1 = 

and there are 20 top-dimensional minors: 



M123 




/11512/113 


— /ll5l3^12 


- 511/12/113 + 511/13/112 + /ill/l25l3 - 


/II1/13512 


M124 




/11512/122 


— /11522/112 


— 511/12/122 + 511/22/^12 + /ill/12522 — 


/II1/22512 


M125 


— 


/11512/I23 


— /11523/I12 


— 511/12/123 + 511/23/^12 + /ill/12523 — 


/II1/23512 


Mi26 




/l 1512/^33 


— /ll533^12 


— 511/12/133 + 511/33/112 + /111/12533 — 


/ll 1/33512 


M134 


= 


/l 1513/^22 


— /11522/I13 


— 511/13/I22 + 511/22/I13 + /I11/13522 — 


/I11/22513 


M135 


= 


/11513/I23 


— /11523/I13 


— 511/13/I23 + 511/23/I13 + /II1/13523 — 


/II1/23513 


Mi36 




/ll5l3^'33 


^ /ll533^13 


^ 511/13/133 + 511/33/113 + /l'll/l3533 ^ 


/''I1/33513 


M145 


= 


/11522/123 


— /l 1523 ^22 


— 5I1/22/I23 + 5I1/23/122 + /1I1/22523 — 


/II1/23522 


Mi46 


= 


/11522/I33 


— /l 1533/122 


— 5I1/22/I33 + 5I1/33/I22 + /II1/22533 — 


/II1/33522 


Mi56 




/l 1523/^33 


~ /l 1533 ^23 


— 5I1/23/133 + 5I1/33/123 + /1I1/23533 — 


/ll 1/33523 


M234 


= 


/l25l3/i22 


— /12522/I13 


— 512/13/I22 + 512/22/I13 + /I12/13522 — 


/I12/22513 


M235 




/12513/I23 


— /l2523/il3 


— 512/13/I23 + 512/23/I13 + /I12/13523 — 


/I12/23513 


M236 




/l25l3'l33 


— /12533/I13 


— 512/13/I33 + 512/33/I13 + /I12/13533 — 


/I12/33513 


M245 




/l2522/i23 


— /l2523^22 


— 512/22/123 + 512/23/122 + /112/22523 — 


/I12/23522 


M246 




/12522/I33 


— /12533/I22 


— 512/22/I33 + 512/33/I22 + /I12/22533 — 


/I12/33522 


M256 




/l2523/i33 


— /l2533/i23 


— 512/23/I33 + 512/33/I23 + /I12/23533 — 


/I12/33523 


M345 




/l3,922/i23 


^ /l3523^22 


^ 513/22/123 + 513/23/^-22 + /113/22523 ~ 


/113/23522 


M346 




/13522/J33 


— /13533/I22 


— 513/22/I33 + 513/33/122 + /I13/22533 — 


/I13/33522 


M356 




/l3523/i33 


— /l3533/i23 


— 513/23/I33 + 513/33/I23 + /I13/23533 — 


/I13/33523 


M456 




/22523/i33 


— /22533^23 


— 522/23/133 + 522/33/123 + /122/23533 — 


/I22/33523 



This time there is again a single Plucker relation 
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-M123M456 - M124M356 + M125M346 - M126M345 + M134M256 

-M135M246 + Mi36M245 + M145M236 - M146M235 + M156M234 = (18) 

In this respect, the case n = 3, r = 2 is analogous to the case n = 2, r = 4. Be it coincidence or not, this is 
not the only parallel between these two cases: see the chapter 5 about integral discriminants. Resultant is a 
homogeneous polynomial of degree 4 in minors: 

R3\2{f,9, h} = M123M145M156M456 - M123M145M236M456 - M123M145M345M356 + M123M145MI46 - 

Mi23M24gM456 + M123 M146M235M456 + M123M146M245M356 - 2M123M146M246M346 - M123M156M234M456 + 

Mi23Mi56M|4g + M123M234M246M356 + M123M234M345M356 - M123M234MI46 - M123M235M245M356 + 
M123M235M246M346 + M123M236M245M346 - Ml23M236M|4e - M123M236M246M345 - M124M126M156M456 + 
M124M126M246M356 + M124M126M345M356 - M124M126MI46 + M124M135M236M456 - M124M135M246M356 + 
M124M136M146M456 - M124M136M235M456 + M124M136M246M346 - Mi24Mf46M356 + M124M146M156M346 + 
M124M146M235M356 - M124M146M236M346 - Mi24Mf56M246 - M124M156M235M346 + M124M156M236M246 + 
Mi24Mi56M236Af345+Ml25Mi26Afi46M456-Ml25Ml26M245A^356+Mi25Mi26M246M346+Mi25Mi34Mi56M 
Mi25Afi34Af236M456 - A/125 A/134 A/345 Af356 + A/125 A/134 A/|4g " A/125 A/135 A/146 A/456 + M125M135M245 A/356 - 
A/125A/135A/246A/346 + A/125A/136M234A/456 - A/125A/136A/245A/346 + A/125A/136A/246A/345 + A/i25M24gM346 " 
A/l25A/l46A/l56A/246 - A/125 A/146 A/156 A/345 - A/125A/146A/234A/356 + A/125A/146A/236A/246 + A/125 A/fsg A/245 + 
A/125A/156A/234A/346 - Afl25 Afl56Af236Af245 - A/f26Afl45A/456 + A/fag A/245 A/346 " Afi26-^246 " A/f26-^246Af345 - 
A/l26 A/134 A/l46A/456+A^126A^134A^235A/456-A^126A^134A^246A^346+A/l26A/l35 A/145 A/456-A/126 A/135 A/234 A/456 + 
Mi26Mi35Af|4g + A/126A/245M356 - 3A/i26Mi45Afi46Af346 + Afi26Af 145 A/156 A/246 + A/126A/145A/156A/345 + 
2A/i26A/i24gA/246 + M126 A/124 g A/345 " 2 A/126 A/146 Ari56A/245 + A/126 A/i46A/234A/346 " A/126 A/146 A/235 A/246 " 
Mi26Ari56A/234A/246 " A/126A/156A/234 A/345 + A/126 A/156 A/235 A/245 - Ml^^Mir,QMir,f^ + Mf^^PhieMs^Q + 
A/f34A/345Ar356 " Mfr^^Ml^^ + A/134A/135 A/146 A/456 " M134M135M245M356 + A/134A/135 A/246A/346 - 
A/134A/136A/145A/456 + Afl34Afl36Af245Af346 - A/134A/136 A/ljg - Afl34A/l36Af246Af345 + 2A/134A/145 A/146 A/356 - 
A/134A/145A/156 A/346 - A/l34A/l45 A/235 A/356 + A/i34A/i45 A/236 A/346 - 2A/i34A/i46-^346 + 3A/i34A/l46A/l56A/246 + 
A/l34A/i46A/i56Ar345 + A/134M146M235 A/346 " A/134 Ari46Ar236Ar246 - Ari34 Ari46 Ar236 Ar345 - Mi34Mi25g A/245 - 

Mis^M^^^Msse + Mi^sMusMueMaie + Mi3.^Mu^M23.iM35Q - Mi35Mu5M23eM2ie - Mi35Mf^QM246 - 

Ari35Ari46Ar234A/346 + A/i35 A/i46 M236 M245 + MiseMf^^M^ie - A/i36Ari45Ari46A/246 - Mis^MusMueMsAS - 
Ari36A/i45A/234A/346 + A/136A/145 A/235 A/246 + Mi36M^4qM245 + A/136 A/146 Ar234Ar246 + A/136A/146A/234 A/345 - 
A/l36A/l46 A/235 A/245 - A/f45-^156 + A/f45^156 A/236 + 2Afi45 Aff46-^156 " A/145 A/f46 A/236 - Afl45 Afl46 Afl56 Af235 + 
A/l45Aff56-^234 - Mf^^ + M^^qM235 - M^^eMi5eM234 + Af256Afl23Afl46Af345 - A/256 A/123A/156 A/245 - 
Ar256Ari23A/234A/346 + A/256A/123 A/236 A/245 - A/256A/124A/126A/346 + A/256A/124A/135A/346 " 
A/256Ari24Ari36Ar345 + A/256A/124A/146A/156 " Ar256 Ari24 Ari46 Ar236 + Ar256Ari26^245 + A/256 A/126 A/134A/345 - 

A/256Ari26Ari35Ar245 " A/256M126 A/145 A/i46 + Ar256 Ari26 Ari46 Ar234 - Ar256 Arf34 Ar346 + A/256A/134A/136A/245 - 

A/256A/134A/145A/156 + A/256A/134A/145A/236 - M256Mi34Ml^Q + M256A/135A/145A/146 - A/256A/136A/145A/234 

This time it is not a determinant of some matrix (or, at least, such matrix is not yet found). Instead, it is a 
PfafEan of the following 8x8 antisymmetric matrix, made of the Plucker coordinates: 
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/ 



-M356 



-^346 *^146 



M"l56 



Ml36 



Ml 26 



A^236 A^i56 — M236 



M234 - M145 



M134 - M125 



Mi26 - M135 



Mi46 - M235 



-"256 -M"236 -M246 -M234 + M145 M126 -^^134 + ^^125 

\ -M346 -M156+M236 -M345 M234 -Ml26 + **135 Ml34 -Mi4e + M235 



(19) 



The PfafRan is equal to 



-R3|2{/)5i^} — pfaffj4 — Cii... 13^1112^1314 ^»5»6^i' 



(20) 



where e is the completely antisymmetric tensor (simply the antisymmctrisation over ii, . . . The Pfafhan 
is also the same, as the square root of determinant, since for antisymmetric matrices determinant is always an 
exact square. Observation (20) was done in [13], where the variables Mijk are being interpreted as coordinates 
on a Grassmanian manifold and then, using advanced homological algebra, this and similar formulas are being 
obtained. Generally speaking, resultant is expected to be a polynomial of degree r"~^ in the Plucker coordinates 
-^ji...i„- It is an open problem to find an explicit formula for this polynomial. 

3.3 Formulas of Sylvester + Bezout type 

By combining the techniques of Sylvester and Bezout, it is possible to sufficiently enlarge the number of cases, 
where resultant is expressed as some polynomial (usually determinant) without any divisions. The basic example 
here is n = 3, r = 2. In this case, the system of quadratic equations has a form 



f{xi,X2, X3) = fnxj + h2XiX2 + /i3a;ia;3 + f22xl + f2sX2Xs + hsx% = 
g{xi,X2,xs) = ffiia;! +5fi2a;ia;2 +5i3a;ia;3 +5f22a;2 +523a;2a;3 +533a;| = 

h{xi,X2,X3) = hnxl + hi2XlX2 + hi^XiXs + h22X2 + h2iX2Xs + /i33a;3 = 



(21) 



It can be represented as a 3 x 6 rectangular matrix. Obviously, another 3x6 matrix is necessary to complete 
a square 6x6 matrix. This complementary matrix is constructed quite elegantly, making use of Jacobian 



J{xi,X2,X2,) = det 



dxi 



, h = f,f2 = 9,h = h 



which is a homogeneous polynomial of degree 3 in xi,X2,X3 and at the same time of degree 3 in coefficients of 
/, g and h. Jacobian can be expressed through Plucker variables Mijk as follows: 
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J{xiX2Xz) = Miizxl + (-2Mi34 + Mi2b)xlx2 + (2Mi26 - Mi35)a;?a;3 + (-M234 + 2Mia-,)xixI + (-M235 + 
4Mi46)a;ia;2a;3 + (-M236 + 2Mi56)a;ia;§ + M245a;i + (2M246 + M345)a;ia;3 + (M256 + 2M346)a;2a;i + Magga;! 

The derivatives of J are homogeneous quadratic polynomials in xi,X2,X3, just like the equations f,g and h: 

^ fll fl2 fl3 f22 /23 /33 ^ 







9 




h 




aj 




aj 




Ui 





v 



3Ml23 -■»Mi34 + 2Mi26 -2M135 + 4M126 - A'/234 + 2M145 _M236 + 4Mi4g -M236 + 2Mi5g 

-2M134 + M125 4Mi45 - 2M234 -M235 + 4Mi45 3M245 2M345 + 4M246 ^256 + 2M34e 

2Mi26 - J^fl35 -M235 +*J^146 -2M236 + 4M156 M345 + 2M246 2M256 + 4M345 3M35e \ 2 / 

The matrix in the right hand side is half Sylvester, half Bezout. Its determinant is equal to the resultant: 



fll fl2 fl3 f22 f23 f33 

911 912 913 922 923 933 

hii hi2 /113 h22 ^23 /133 

3M123 -4^134 + 2M125 -2Mi35+4Mi26 -^234 + 2M145 -M235 + 4Mi4g -M236 + 2Mi5e 

-2M134+M125 4Ml45 - 2M234 -M235 +4Mi4e 3M245 2M345 + 4M246 M256 + 2M346 

2Mi26 - M135 -M235 +4Mi4e -2M236 + 4Mi5e M345 + 2M246 2M256 + 4M345 3M35g 



R3\2{f,g,h} = 



(22) 



In a compact notation this formula can be written as 



R3\2{f,g,h} = det {f,g,h 



dJ_ dJ_ dJ 

dxi ' dx2 ' 0x3 



(23) 



where det{. . .} stands for the determinant of the matrix, obtained by expanding the polynomials inside brackets 
in the monomial basis. The other possibility to design a square matrix is to take 9 polynomials Xifj, which are 
homogeneous polynomials of degree 3, and complement them with the Jacobian itself, not its derivatives. In 
this case one gets a 10 x 10 matrix in the monomial basis: 



R3\2{f, 9, h} = det {xif, xig, Xih, X2U X2g, X2h, x^f, x^g, x^h, J} 



10x10 



(24) 



According to this method, it is necessary to make the matrix square by adding a number of certain polynomials 
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to the Sylvester block. Of course, the most non-trivial part of the method is construction of these polynomials. 
For higher r > 2, these polynomials are constructed with the help of 



J{xi, . . .,Xn,Pl, ■ . . ,Pn) = det ( + ^ J Pk 



dxj r — ldxjdxk {r — l){r — 2) dxjdxkdxi 



PkPi + ■ ■ 



Jacobian-like quantity J' can be viewed as a series in the additional variables pi, . . . ,p„: 

J7^(iri , . . . , X^i^Pl 5 ... 5 Pn) — J {x\-, ^ • • 1 Xji^ -\- Jiixi^ . . . ^ XYi)Pi + "^^Jij^X-i^ . . . , X^^PiPj + . . . 

Note, that J{xi,...,Xn) is nothing but the ordinary Jacobian and Ji{xi, . . . , x„) are its derivatives, while 
quantities Jij{xi, . . . , Xn) are ahcady now and non-trivial (in particular, they arc not equal to the second 
derivatives of J). In the next case n = '6,r — S the Sylvester-Bezout formulas take form 



-R3|3{/i, /2, /a} = det {xiXjfk, Ji} 



(25) 



and 



-Rsiai/i, /2, /s} = det {xifj, Ju} 



15x15 



(26) 



Just like in the previous case, there exist two different formulas. In the next case n = 3, r = 4 we have 



^3|4{/l, /2, h} = det {XiXjXkfu Jpq} 
obXoD 



and 



(27) 



^3|4{/l, /2, h} = det {XiXjfk, Jlpq} 



28x28 



(28) 



Generally, for n = 3 and r = any, we have two formulas of Sylvester-Bezout type: 





-R3|r{/l) /2, /a} 


= det jxj, . 

r(2r-+l)xr(2r+l) ^ 


• • Xi^_^fj, Jki. 


..fer-2} 


and 












-R3|r{/l)/2,/3} 


= dct \xi. . 

r(2r-l)xr(2r-l) 


■ ■ Xi^_2fj, Jki. 


..fer-l } 



(29) 



(30) 



As for higher dimensions n > 3, an explicit formula is known only for n = 4, r = 2 where it has a form 



-R4|2{/i, /2, /a, /4} = det {xifj, Jk} 



20x20 



(31) 



For higher n and r, such formulas are still an open direction of research, see [11] for more details. 
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3.4 Formulas of Schur type 

Resultants also admit a very different kind of formulas, related to certain analytic structure of the logarithm of 
resultant, described in [12]. By an analytic structure we mean that logarithm of resultant of a system 

fl{xi,...,Xn) = 
f2{xi,...,Xn) = 



depends on polynomials fi, - ■ ■ ,fn in a very specific way, which is much simpler to understand than for the 
resultant itself. Two formulas are used to express this statement. 
The first is the operator formula 





OO 




OO 






log i?,^,...,,„(7-/) = 








• log det(/ - A) 


A=0 



(32) 



where A is an auxiliary n x n matrix, and /i, /2, • • • , /n are differential operators 

/i = /i 

/2 = /2 



d 


d 


d 




dAi2'' 




d 


d 


d 


dA2i 


dA22 ' " 


■ " ' dA2n 



(33) 



fn — fn 



( d d 



d 



\dAnl' dAn2'"'' OA 

associated with polynomials /i, • • • , /n- The second is the contour integral formula: 




(34) 

where the contour integral is {n — l)-fold and functions fi are taken with arguments fi{l, Zi, . . . , Zn-i)- The 
contour of integration is assumed to encircle all the common roots of the system 

/fc(l,2;i,...,2;„_i) = 0, k = 2,...,n (35) 

see [12] for more details. An explicit formula for logi? follows after several iterations of (34): for n = 2 

logi?,,,, {/i,/2} = ^log/i(l,2)(ilog/2(l,z) + rilog/2(0,l) (36) 
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for n = 3 



logi^rirsrs {flJlJs} = j j log /i (1, Zi, ^2)^ log /2(1, ^1 , ^2) A d log /g (1, Zi , Z2) + 

n j log /2(0, 1, z)d log /3(0, l,z) + riri log /3(0, 0, 1) 



and so on. In this way, logarithm of the resultant is expressed through the simple contour integrals. To obtain 
R from log R, we exponentiate the latter: for example, for n = 2 



-Rr-ira {fl,h} = exp (^j log /i ( 1, z)d log /2 (1 , z) + Ti log/2(0, i; 



(37) 



Using this formula, it is straightforward to express the resultant of two polynomials through their roots: if 



h{l,z) = aY[{z-ai) (38) 



f2{l,z) = hY[{z-pi) (39) 

i=l 



then 



Y.Y.f los^-^ - "i)'^log(^ - l^o) + ri loga + n log6 = afW ViX\!y^i ' (^0) 

Exponential formula (37) and its direct analogues for n > 2 have several advantages. The most important 
are simplicity and explicitness: such formulas have a transparent structure, which is easy to understand and 
memorize. Another one is universality: they can be written for any n, not just for low values. However, there 
is one serious drawback as well. Resultant is known to be a homogeneous polynomial of degree 

di = 

in coefficients of the i-th equation /j, and of total degree 

d = di + ... + dn = nri ...r„f — + — + ... + — 



n ri rr, 



in coefiicients of all equations, but in exponential formula the polynomial nature of the resultant is not explicit. 
Therefore, it is not fully convenient for practical calculations of resultants (although it is conceptually adequate 
and can be used in theoretical considerations). 
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To turn these formulas into a computational tool, one can apply (34) to a shifted system 



f fiix) = ixiY' - Xihix) 

hix) = ix2Y' - X2f2ix) 
. fn{x) = {XnY" — Xnfnix) 



(41) 



and then expand the shifted logarithms in the integrands into Taylor series in powers of the "spectral parameters" 
Aj. Then, one obtains the following series expansion 



logi?{/i,...,/„} 



oo oo 

E--E 

fci=0 A:„=0 



Tkik2...k„{fl, ■ • ■ , fn} ■ ^1^X2^ ■■ - 



(42) 



(a minus sign is just a convention) where particular Taylor components Tk-^k2...kn homogeneous polynomial 
expressions of degree ki in coefficients of fi. We call them traces of a non-linear system fi{x),. . . , fn{x). By 
exponentiating, one obtains the shifted resultant: 

(00 oc \ 

fci=0 fc„=0 / 



The original resultant . . . , fn} is equal to (— l)** times the coefBcient of Xf^X2^ . . . A^" in R{fi, ■ ■■ , fn}- 

We can extract it, expanding the right hand side in powers of A,: 



exp 



fei=0 



. . . x!i" 



fen=0 



E- 

ki=0 



■ E '^klk2 
kn=0 



•A^A2^ 



• A„ 



Such power scries expansion of exp {S{x)), where S{x) is itself a power series (perhaps, of many variables), is 
often called a Schur expansion. It implies the following relation between Taylor components: 



ki + ... + k„ 

E 

m=l 



(-ir 



m! 



E 



^ V2 



V1+V2 + - ■■+Vm=k 



where the sum is taken over all ordered partitions of a a vector k = (fci, /c2, • • • , fcn) into m vectors, denoted 
as Vi, . . . ,Vm. An ordered partition is a way of writing a vector with integer components as a sum of vectors 
with integer components, where the order of the items is significant. Polynomials Vkik2...k„{T) are often called 
multi-Schur polynomials. Several first multi-Schur polynomials are given by 
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^1,0 — ^2^1,0 
7^2,0 = -T2,o + TlJ2 
V2,i = -T2,i + T2,oTo,i + Ti,oTi,i - (43) 
'^2,1,0 — ^72,1,0 + 72,0,0^0,1,0 + 2^1,0,02^1,1,0 ^ 2\^,o,o2o,i,o/2 
"Pi, 1,1 = — ri.i.i + I'l.o.oTo.i.i + 2o,i,oT'i,o,i + I'l.i.olo.o.i — 2o,i,oIo,o,iT'i,o,o 

Resultant is a certain multi-Schur polynomial of traces: 

R{fl, ...,fn} = i-lf rd,,...,dn = E j,, E • • • (44) 

'""^ iri+jr2+-..+i'Tn=d 
The explicit formula for traces, obtained in [12], has a form 

1 " 



klk2...kn ^ 2<i,j<n . , 

Tij 1=1 



where 



{f)ji,j2,-,3„ = coefBcient of a;f a;^^ • • • a;^" in a;2, . . . , a;„)'= 



and the sum goes over all non-negative integer n x n matrices Vij such that the sum of entries in any i-th row 
or i-th coloumn is fixed and equals rj/cj. Formula (45) is valid for positive fci, . . . , A;„. The cases when some 
ki = 0, also make no difiiculty. If some fcj = 0, then 



Ifei,fe2,...,fe„{/l) /2, • • • , fn} — ri Tki,k2,..;ki-i,ki+i,--;kn{h^ f2, • • • > /i-1) /i+1) • • • ) fn}\^.^Q (46) 



where the trace in the right hand side is taken in variables xi, . . . , Xi-i, x^+i, . . . , a;„. It is an easy exercise to 
prove this statement, making a series expansion in (34). Similarly, operator approach gives [12] 



T If f fX- ^l(/l)"^ Mfn)"" trA-^k,+...+r^k„ 



(47) 

A=0 



Together, formulas 44 and 45 (or, alternatively, 47) give an explicit algorithm to calculate resultants without 
any polynomial divisions. The absence of polynomial divisions should increase computational speed of the 
algorithm, if compared to algorithms with division such as Koszul complex. 
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To illustrate these general formulas, let us consider a particular calculation for n = 3 and r = 2, where the 
system of equations again has the form 21. The traces are calculated with (47), using differential operators 



/ = /ll(^ +/22h:^ +/23 ./.. +/33^ ^ 



5^11; ■'"'5^119^12 ■''"5An5^i3 ■'""V5^i2y ' 9^129^13 V^^is 



^ = ^"'91^^ +^^^9A2i9A22 +^"9A2i9A23 + ^Sl^J + ^4^^ + 9^23 

" = "11 TTl — + "-12 771 ?n 1- "13 771 ?n 1" "22 TTl — + "23 771 ^71 1" "33 



^9^31/ 9^319^32 9^319^33 V^^32/ 9^329^33 V^^33 

and a few first traces are 
T\ = 4/ii + 4fif22 + 4/133 

Tl = ^flx + 4512/12 + 4/I13/13 + 8311/22 + 8/I11/33 + 4fi(|2 + 8333/122 + 4fif23/i23 + 4/I33 

^3 = 4/i\ + 6512/11/12 + 6/113/11/13 + 12511/11/22 + 12/111/11/33 + 6511/1^2 + 6/111/1^3 + 3523/113/12 + 
6533/112/12 + 6512522/12 + 3513/123/12 + 3523/112/13 + 6/113/133/13 + 6513/122/13 + 3512/123/13 + 6512/22 + 
12511522/22 + 6513/113/22 + 12533/111/22 + 6523/111/23 + 3512/113/23 + 9513/112/23 + 6511/123/23 + 6/113/33 + 
6512/112/33 + 12511/122/33 + 12/111/133/33 + 12522533/122 + 45I2 + 6523/123/133 + 6522523/123 + 
6523/122 + 6533/1^3 + 12533/122/133 + 4/133 

and so on. The resultant has degree 12, and is expressed through traces as 



i?3|2{/,5,/l} = Pl2{^} 



"1^^^^ + + ^^^^^ + + i^^^^ + + TT^^^" " 3^^' " ^^^^^^^ " ik^^'^^ ~ 

T2T4T6 T2T3T7 TlTg TiTsTe - —T1T4T-1 - — TiTaTg T1T2T9 T?Tiq + 

100 ' 48 42 64 ^ 30 28 24 18 20 ^ 

+ ^,T,TiT, + ±-^Tin + ^,Tinn + ^^Tin + ^^T^nn + ^TiTIT, + l^T\T,T,n + 

36 56 ^ 100 ^ ^ 48 ^ 42 ^ 32 ^ 54 ^ 864 ^ ^ 

^ 1536^2 ^ 324^^^^^^ ^ 96^^^2 T3T4 - — T1T2 T5 - — T.T, - J^T.T^T, - -T, - 

96^1 ^2 ~ 120^1 " 108^1 ^^^'^ " 84^1 ™ ~ 192^' + 46080^^ + TT52^^^2 ^3 + ^T, T, T, + 
384 1 2 4 T 1 3 144 1 ^ 240 ^ ^ 768 ^ ^ 360 ^ ^ ^ 288 ^ ^ 840 ^ 

7680^1 ^2 " 864^1 ^2 " 864 ^1 768^' ~ Tiio^^ ~ 1200^1 ^'^^ 4320^^ + 

-^TfT^ + ^— T|r3 + Ti^Tl + ^— Tf T2T4 + —^TlTr, Tf T| 5— r7T2T3 - 

9216 ^ ^ 2880 1 2 3 T ^2ggo 1 3 ^ g^gg 1 ^ 4 -r 25200 ^ 34560 ^ ^ 30240 ^ 

T^Ta + -^^ T-i8ji2 1 ^ T^Ti ^ T^^Ti A 7-'i2 

161280 ^ 322560 ^ ^ 1088640 ^ 7257600 ^ 479001600 ^ 

This is a rather compact formula for -R3|2- By substituting the expressions for traces, we obtain a polynomial 

expression in /, 5 and h, which contains 21894 monomials and takes several dozens of A4 pages. 
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Figure 1: Determinant of matrix Sij, represented as a diagram of tensor contraction. Black 2-valent vertices 
represent tensor S, white n-valent vertices represent tensor e. 

4 Calculation of discriminant 

Discriminant is a particular case of resultant, associated with gradient systems of equations. Since discriminant 

is a function of a single homogeneous polynomial (instead of n homogeneous polynomials in the resultant case) 
it is slightly simpler than resultant, has bigger symmetry and some specific methods of calculation. These 
discriminant-specific methods are the topic of this section. 

4.1 Discriminant through invariants 

Let S{xi, . . . ,Xn) be a homogeneous polynomial of degree r in n variables. Discriminant of 5 is a polynomial 
SL{n) invariant of S. As such, it should be expressible as a polynomial in some simpler basis invariants, given 
by elegant diagram technique suggested in [1]. The simplest example is of course the determinant, given by a 
single diagram at Fig.l. The simplest non-trivial (i.e. non-Gaussian) example is a 3-form in 2 variables: 



By dimension counting, there is only one elementary invariant 1^ in this case, given by a diagram at Fig. 2. 
The subscript "4" stands for the degree of this invariant. In this paper we find it convenient to denote the 
elementary invariants of degree k as Ik- It is straightforward to write the algebraic expression for the diagram: 



^(xi, a;2) = Siiixf + 35'ii2a;?a;2 + BSi22Xixl + S222XI 



'ili2i3'-'jlj2j3 



SkMs S1M3 ^'"'^ e'^'^ e'^'^ e'^'^ 



Evaluating this sum, one gets the following explicit formula for 



I4 — 2811182'^ 



125'iii<S'ii2<S'l22'S'222 + 8<S'iii<S']^22 + ^Sii2^222 — 6>S';^]^2'S'l22 



which is nothing but the algebraic discriminant £>2|3 of S: 




(48) 
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Figure 2: The degree 4 invariant I4 of a 3-forni in 2 variables, represented as a diagram of tensor contraction. 
Black 3-valent vertices represent tensor S, white 2-valent vertices represent tensor e. 



The next-to-simplest example is a 4-form in 2 variables, which can be written as 

S{xi,X2) = S'liiiXi + 4Siii2xlx2 + 6Sii22xlxl + 4S'i222a;ia;2 + 5'2222a:;2 

By dimension counting, there arc two elementary invariants in this case. They have relatively low degrees 2 
and 3, denoted as I2 and and given by diagrams at Fig. 3 and Fig. 4, respectively. Looking at the diagrams, 
it is straightforward to write algebraic expressions for h^h'- 




Figure 3: The degree 2 invariant I2 of a 4-form in 2 variables, represented as a diagram of tensor contraction. 
Black 4-valent vertices represent tensor S, white 2-valent vertices represent tensor e. 



— q. . . . Q. . . . ,«2J2,*3J3,UJ4 
J-2 — '^JlJ2»3»4'^Jlj2j3j4'^ e t t 

_ C. . . . C. . . . C, , , , Jljlj2j2j3kljik2 j3k3 j4k4, 
-'3 — '^JlJ2»3»4'^Jlj2j3j4'^'ClK2'C3K4'^ t t t t t 
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Figure 4: The degree 3 invariant I3 of a 4-forni in 2 variables, represented as a diagram of tensor contraction. 
Black 4-valent vertices represent tensor S, white 2-valent vertices represent tensor e. 



Evaluating these sums, one gets the following explicit formulas for hjh- 

I2 = 25111152222 — 85iii25i222 + 6^1122 
I3 = 65llll5ll2252222 — 65llll5i222 — 651112^2222 + 125lll25ll225l222 — 65ii22 

The algebraic discriminant D2\4,, just like any other 5L(2)-invariant function of 5, is a function of 12, Iz'- 

8'5'llll52222 ~ 965]^]^]^]^5iii25i22252222 ~ 1445fiii5f 

122'5'2222 + 43 25]^]^]^]^5ii225]^222'5'2222 



i?2|4 = II - Ml 



— 2165iiii5i222 + 4325iiii5iii25ii225|222 ~ 485iiii5iii25i22252222 — 14405iiii5iii25ii225i22252222 
+ 6485iiii5fi2252222 + 8645iiii5iii25ii225i222 + 8645iii25ii225i22252222 + 2885iii25ii225i222 

— 4325iiii5ii225i222 — 2165fii25|222 ~ 5125iii25i222 — 4325iii25ii22 52222 (49) 

Our third example is a 5-form in 2 variables, which can be written as 

5(a;i,a;2) = 5iiiiia;f + 55iiii2a;ia;2 + 105iii22a;?a;2 + 105ii222a;?a;2 + 55i2222a;ia;| + 5222222:2 

By dimension counting, there are three elementary invariants in this case. They have degrees 4, 8 and 12, 
denoted as I4, /g and /12 and given by diagrams at Fig. 5, Fig. 6 and Fig. 7. 
Looking at the diagrams, it is straightforward to write an expression for I4 

T, — q C <?,,,,, <?,,,,, filhfi2j2j3j3fiiki i^kz jili kill k2l2fk3l3 

-14 — '-'»l»2»3»4»6'-^Jlj2j3j4j6'-^'Sl'C2'C3«4K5'-^(l'2'3l4'5'^ CCCCCCCCC 



and equally straightforward to write expressions for I%,Ii2- Evaluating the contraction, one gets a formula 
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Figure 5: The degree 4 invariant I4, of a 5-form in 2 variables, represented as a diagram of tensor contraction. 
Black 5-valent vertices represent tensor S, white 2-valent vertices represent tensor e. 




Figure 6: The degree 8 invariant 1% of a 5-form in 2 variables, represented as a diagram of tensor contraction. 
Black 5-valent vertices represent tensor S, white 2-valent vertices represent tensor e. 



^'S'l 11 11 522222 ~ 20S'iiiii6'iiii2<S'i2222'S'22222 + 8S'iiiii6'iii225'ii2225'22222 + 
325iiiiiS'iii22'S'l2222 ~ 24S'iiiiiS'ii222'S'l2222 + 32S'im2'^11222'S'22222 + 18'S'im2'^12222 ~ 
24S'iiii2iS'm22'S'22222 — 1525iiii2S'iii225'ii222'S'l2222 + 96S'iiii2S'ii222 + 965^1122^ 12222 — 645^1122^ 11222 



and similar formulas for I8,li2 - they are quite lengthy and we do not present them here. The algebraic 
discriminant -D2|5) just like any other S'L(2)-invariant function of S, is a function of I4,ls,li2' 



£>2|5 = II - 64/8 



(50) 



Our final example is a 3-form in 3 variables, which can be written as 

5(a;i,a;2,a;3) = 

Suixl + 3Su2xlx2 + SSii3xlx3 + S222X2+3Si22Xlxl+3S223X2X3 + S333X3 + SSi33Xixl + 3S233X2xl + 6Si23XlX2X3 

By dimension counting, there are two elementary invariants in this case. They have degrees 4 and 6, denoted 
as I4, Iq and given by diagrams at Fig. 8, Fig. 9. 
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Figure 7: The degree 12 invariant /12 of a 5-form in 2 variables, represented as a diagram of tensor contraction. 
Black 5-valent vertices represent tensor S, white 2-valent vertices represent tensor e. 

Looking at the diagrams, it is straightforward to write the algebraic expressions for I4 and Iq: 

T. — Q. . . q . . . c, , , q, , , Mhki J232I2 Mksh Mk2h 

T„ — q. . . q . . . q, , , q, , , q q Jikih i2j2S2jik2mi I2m2k3 msssjs Jsissi 

-'6 — '-'lll2i3^Jlj2j3'~'klk2k3'-'lll2l3^mim2m3'-'SlS2S3'^ c c c c c 

Evaluating these sums, one gets the following explicit formulas for hjh- 

I4 = 63^23 — 125'i22'S'i23'S'i33 + 6<S'^22'5'i33 + 6S'ii35i23>S'l335'222 — 125'ii35j23'S'223 — 65ii35i22'S'l33'S'223 + 
185ii35i22<S'l23'5233 — 6S'ii35i22'^333 + 6<S'ii3<S'|23 — 6S'^i3S'222'S'233 ~ 6S'ii2<S'i33<S'222 + 18<S'ii25'l23'S'l33'5223 ~ 
12S'ii2S'i23<S'233 — 6S'ii2<S'l22<S'l33'S'233 + 6S'ii2S'l22S'i23'S'333 — 66'ii2iS'ii36'223'S'233 + 6S'ii2iS'ii3S'222'S'333 + 
66'^i2'^233 ~ 6S'ii2'S'2235'333 — 6iS'iiiS'i335|23 + 65'iii5'i:i3S'222'5'2;i.'i + 6S'iiiS'i23S'223'S'233 — 66'iiiS'i23S'222'S'333 — 

6S111S122S233 + 6S'iiiS'i22'S'223'S'333 

le = 485123 ~ 144S'l22'S'i23S'l33 + 144S'i22'S'l23'S'l33 — 48S'i22'^133 + 725ii3S'i235'l336'222 — 144S'ii3S'f23<S'223 — 
725ii3S'i22<S'l23'S'i33<S'222 + 725ii3S'i22<S'i23'S'l33'S'223 + 216S'ii3S'i22<S'i23'S'233 + 72S'ii35i22'S'l33'S'223 — 

216S'ii3S'i22'S'l23'S'l33'S'233 " 72S'ii3S'j22'5'i23'S'333 + 725ii3S'f22'S'l335'333 + 185'ii3S'i33 5*122 ~ 
'''2'S'^]^3S'l23'5'i335222'S'223 + 1445^j^3S'^235'|23 ^ 725f]^35^23 5*222 '5'233 + 72S'^]^35i22'S'i33S'|23 " 
36'5ll3'5'l22'S'l33'S'222<S'233 — 2165^j^3S'i22S'l235'223'S'233 + 144S'J]^35i22'5'l23'5'222<S'333 + 162S'^]^35^22 "5*233 ^ 
1445j]^3S'^22'S'2235'333 — 485f]^3<S'223 + 72S'f]^35'222S'223'S'233 ~ 24S'j']^35|22'S'333 ~ 72S'ii2S'i23'S'j33S'222 + 
2165ii2S'i23'S'i335223 — 14451 12'S'l23'S'233 + 72S'ii2<S'i22S'i33S'222 — 2165i 12'S'i225'i23'S'i335223 + 
125i225;i23'^1335233 + 725112512255*235333 + 725i 12 5^22 5133 5233 " 725ii25j225l235i335333 — 
365ii25ii35i3352225223 — 2165ii25ii35i235i335|23 + 3605ii25ii35i235i3352225233 + 725ii25ii35i23 5223 5233 — 
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Figure 8: The degree 4 invariant of a 3-forni in 3 variables, represented as a diagram of tensor contraction. 
Black 3-valent vertices represent tensor S, white 3-valent vertices represent tensor e. 



2165ii25i 13'5'i23'S'222'S'333 + 365ii2S'ii3S'i22'S'l335223'S'233 — 108S'ii2S'ii3S'i22'5l33S'222'S'333 — 
2165ii2<S'ii3<S'i225'l23'S'233 + 3605ii2S'll35'i22'S'l23'S'223'S'333 — 365ii2S'll3S'j22'S'233*S'333 + 72S'ii25'f]^3S'|23<S'233 
144S'ii2S'ii3S'2225'233 + 72S'ii2S'ii3S'222'S'2235'333 + 162S'ii2'^133'5223 ~ 144S'ii2'S'i33S'2225'233 — 
216S'fi2'S'l23'S'l33S'223'S'233 + 1446'Jj2'^123'S'l33'S'222<S'333 + 1446'i;^2'5'l23'^233 ~ 72S'^j2'S'l23'S'223'S'333 + 

'''2'5'ii25'l22'5'l335'233 ~ 36S'J]^2'^122'S'l335'223'5'333 " 72S'^j^2'5'l22'S'l23'5'233'S'333 + 185^j2'5'l22'5'333 + 



72 1 2 'S'l 1 3 'S'2 2 3 S'l 



144S'f]^2'5'll35'|235'333 + 72S'f]^2'5'll35'2225'233'S'333 ^ 485f]^2 "5*233 + 725*^ ]^2'S'2235'233'S'333 — 



24<S'ii2'S'222<5'333 — 24511151335122 + 1445iii5i235i3352225223 " 725iii5i235i335|23 " 
2165iii5f235i3352225233 + 725iii5f235223'S'233 + 1205iii5f2352225333 — 1445iii5i225^335|23 + 
725i ll5i225i33 5222S'233 + 3605iii5i225i235i33 52235'233 — 725iii5i225i235'i33 5222'S'333 — 725iii5i22S'i23 5233 — 
2165iii5i22S'i23 5223'S'333 — 1445iii5i225i335233 + 725iii5i225i33 5223'S'333 + 1445iii5i225i23 5233 5333 — 
245iii5i22 5|33 + 725iii5ii35i335223 — 1085iii5ii35i3352225223'5'233 + 365iii5i 13 5i33 5|22 "5*333 ^ 
725iii5ii35i23 5|235233 + 1445iii5ii35i23 5222 5|33 — 725iii5ii35i2352225223'S'333 ^ 365iii5ii35i225223"5'|33 + 
725iii5ii35i225|235333 — 365iii5ii35i22 5222 5233 5333 — 36 5iii5ii25i335|235233 + 725iii5ii25i3352225|33 — 

365iii5ii25i33 5222 5223<S'333 — 725iii5ii25i235223'S'233 + 1445iii5ii25i235|23 5333 — 
725iii5ii25i23 52225'233'S'333 + 725iii5ii25i225233 — 1085iii5ii25i22 52235'2335'333 + 365iii5ii25i22 5222'S'333 + 
185iii5|235|33 — 245iii52235333 — 245iii52225'233 + 365iii52225'2235'2335'333 — 65iii5|22'5'333 

The algebraic discriminant -D313, just like any other 5i(3)-invariant function of S, is a function of I4 and Iq: 



£»3|3 = 327| + 37| 



(51) 



When expanded, discriminant £>3|3 contains 2040 monomials. See the appendix of the book version of [1], where 
it is written explicitly. Formula (51) is a remarkably concise expression of this disriminant through a pair of 
invariants, given by beautiful diagrams Fig. 13 and Fig. 14. It is an open problem to generalize this type of 
formulas - eqs. (48), (49), (50) and (51) - to higher n and r. For more information on invariant theory, see [14]. 
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Figure 9: The degree 6 invariant Iq of a 3-forni in 3 variables, represented as a diagram of tensor contraction. 
Black 3-valent vertices represent tensor S, white 3-valent vertices represent tensor e. 



4.2 Discriminant through Newton polytopes 



A homogeneous polynomial of degree r in n variables 



S{xi,X2, ...,Xn) = 



'ai,a2---a„'''l 




.a-n 
n 



can be visualized as a collection of points with coordinates (ai, . . . , a„) in the n-dimensional Euclidean space. 
The convex hull of these points is a n-dimensional convex polytope, called Newton polytope of S. It is known, 
that different triangulations and cellular decompositions of this polytope arc in correspondence with different 
monomials in the discriminant of S. We do not go into details on this topic, see [3] for more details. 

4.3 The symmetric case 

An important special case is the case of symmetric polynomials <S'(xi, . . . , Xn), i-e, those which are invariant 
under permutations of variables Xi, . . . ,Xn- It turns out that exactly for such polynomials discriminant greatly 
simplifies and factorizes into many smaller parts. A few examples for low degrees of S can be given to illustrate 
this phenomenon. The simplest example is the case of degree two, because it is completely treatable by methods 
of linear algebra. A homogeneous symmetric polynomial of degree two should have a form 



where C2, Cn are two arbitrary parameters and pk = J2i ^i- calculate the discriminant, we take derivatives: 



S{xi, . . . ,Xn) = C2P2 + CupI 



(52) 
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ds 

dxi 



= 2C2Xi + 2Ciipi = 



(53) 



The simplest option is just to write a matrix and calculate its determinant by usual rules: 



2?„|2(C2,Cii) = det 



/ C2 + Cii Cii 
Cii C2 + Cii 



c 



11 



... C2 + Cn y 



Cr'(C2+nCn) 



(54) 



As one can see, the discriminant is highly reducible: factorises into simple elementary constitutients. In fact, 
this is a general property of symmetric polynomials. The next-to-simplest case is a homogeneous symmetric 
polynomial of degree 3 in n variables, which has a form 



S{xi, ...,Xn)= CzPz + C21P2PI + Ciiip\ 



(55) 



and contains three parameters C3, C21, Cm. As shown in [15], the discriminant of this polynomial equals 



(n-1)! 



2?n|3(C3,C721,C7lll) = (B3) 



(n-3)2"-l 



n 



n-2k 
9n 



2 , 4A;(n - k) k\{n - \ - k) 



27n2 



(56) 



where 



Bi = n^Ciii +nC2i +C3, 
S2 =nC2i +3C3, 
B3 = C3, 



(57) 



This formula is remarkably concise and, most importantly, it contains n just as a parameter and is valid in any 
dimension n. This property makes it possible to study the large n asymptotics and various continous limits. 
Moreover, the above examples for r = 2 and r = 3 can be generalized to arbitrary degree. Generally, a symmetric 
polynomial homogeneous of degree r in n variables has a form 



S'(a;i, . . . ,a;„) = ^ CyPy 

\Y\=r 



(58) 



where the sum goes over Young diagrams (partitions) F : fci > . . . > of fixed sum \Y\ = ki + k2 + . ■ ■ = r, 
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and pY = YliPYi- The number of free parameters is equal to the number of partitions of degree r, which we 
denote as P{r). For example, P(2) = 2 corresponds to partitions (2) and (1,1). Similarly P(3) = 3 corresponds 
to partitions (3), (2, 1) and (1, 1, 1). Several first values of P{r) are 



P(r) 



123456 7 8 9 10 
1 2 3 5 7 11 15 22 30 42 



(59) 



The generating function for P{r) is given by 



oo oo ^ 



(60) 



A symmetric homogeneous polynomial of degree r has P{r) independent coefficients Cy . Its discriminant is a 
function of these coefficients, which we denote V^ir (C). According to [15], this discriminant is equal to 





#m! 


{Ml + . 


. + Mr-iy. 


■Dn|r (C) = const • Cr"'" Y[ 


(dM {C))(.r-1) 


! Ml!. 


.Mr-l\ 


Mi + ...-\-Mr_l=n 







(61) 



where the product is taken over all decompositions of n into r — 1 nonnegative parts and #m is the number of 
zeroes among Mi, . . . ,Mr-i. The degree f3„\r is fixed by the total degree of the discriminant 



and 



with 



0„\r = degO^i^ 



Mi + ...+Mr-i=n 



#M\dcgdM (.\h + ... + .\Ir-iy. 

(r-1)! Mi\...Mr-il 



(62) 



I i i<j i<j<k 



(63) 



1 p; 



(i) 



det 

2x2 



1 p; 



U) 




( 1 



det 

3x3 



«) \ 
M \ 



1 yj Pi 

\_i_y, P« / 
/ 1 



det 

3x3 



Vi % 



2 X ' 



1 Vi Vj 
V 1 Vk vl J 



(64) 
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and so on, generally 



( 1 Vil 

1 2/i2 



det 

kxk 



V 



dct 

kxk 



1 Vik 



I 1 2/n 
1 



fc-2 p(n) \ 



P, 



(ik) 



(65) 



Note, that each of quantities P^^ "'''^ is a polynomial of degree r — in variables j/j. They are defined only for 
pairwise distinct upper indices ii,. . . ,ik and are symmetric in them. For k > r they are undefined: one can 
not choose r distinct items out of (r — 1). As one can see, each particular factor dm (S) can be computed as a 
resultant in no more than in (r — 1) variables. This allows to calculate discriminants of symmetric polynomials 
on the practical level even for n >> r. 

5 Integral discriminants and applications 
5.1 Non-Gaussian integrals and Weird identities 

Resultants and discriminants from the previous sections are purely algebraic (or combinatorial) objects - i.e, 
they are polynomials, possibly lengthy and complicated, but still finite linear combinations of monomials. The 
next level of complexity in non- linear algebra is occupied by partition functions, which are usually not polynomials 
and not even elementary functions. Still, they posess a number of remarkable properties and are related closely 
to resultants and discriminants. The simplest partition functions of non-linear algebra are integral discriminants 



-^n|r(<S') = J 



dxi...dxn e-*(^i'-'^") 



(66) 



For quadratic forms, diagonalisation and other linear- algebra methods can be used to obtain a simple answer 
>^n|2('S') = (dots) but for higher degrees r such methods do not work. Instead, to calculate integral 
discriminants one needs to use a powerful technique of Ward identities, which are linear differential equations 



d 



^dSai,...,an dSbi,...,bn dSc 



Jn\r{S)=0, ai + bi=Ci+di 



(67) 



satisfied by integral discriminants. Validity of these identities follows simply from the fact, that differential 
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operator in the left hand side of (67) annihilates the integrand in (66). Making use of the Ward identities, in 
[5] integral discriminants were found in the simplest low- dimensional cases: 



Invariants 



Discriminant D, 



First branch of the integral discriminant J„|^ 



h, I 6 



It - 64/8 



32/| + 3/| 



-1/2 



-1/6 



-1/4 



1 (l/12)i(5/12)i /6/: 



i^o «! (l/2)i 



-1/10 



-1/2 



J_ (3/10),-^J^(l/10)2,+3J(l/10)j /le/sV / 128Ji2 
ij^oi'-f-' (2/5)i+2i(3/5).+2i \IIJ V 3/| 



1/4 - 1 (l/12),(5/12)i 



(1/2), 



327|; 



where 



(a)fc= ^^^-^ =a(a+l)---(a + fc-l) 



Already by definition the integral (66) is S'L(n)-invariant and, consequently, depends only on the elementary 
S'L(n)-invariants Ik, which were described in section 4.1. above. It is necessary to note, that J„\r is a multi- valued 
function: for one and the same n\r there can be several branches, which satisfy one and the same equations 
(67), but have different asymptotics. For example, for n\r = 2|4 integral discriminant has exactly two branches: 



•/i:i(5)=/2-^/SFi( 1/12,5/12 



Gil 
II 



1/4 ^ 1 (l/12),(5/12)i fell 



(l/2)i 



(68) 



t(^)(q\ - T 7~V4 p / 7/12,11/12 

•^214 ('=')- •'3-'2 2i'l I 



-'2 



hi 



3^2 



-7/4 



E 



1 (7/12),(ll/12), /6/| 



(3/2). 



(69) 



For the sake of brevity, only one of branches is presented in the table above - the one with normalisation 
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factor depending only on the invariant of lowest degree 7min. Other branches are omitted. As one can see, non- 
Gaussian integrals posess a nice structure: they are all hypergeometric, i.e, coefficients in their series expansions 
are ratios of T-functions. As shown in [5], in all these cases integral discriminant J„\r{S) has singularities in 

the zero locus of the ordinary algebraic discriminant Dn\r{S). This interesting relation between purely algebraic 
objects (discriminants) and non-Gaussian integrals deserves to be further investigated. 

5.2 Action-independence 

Actually, an even stronger statement is valid: not only 



d d d 



,9Sai,...,a„ dShi,...,hn 9Sci,...,c„ dSdi,...,dn 



da;i . . .da;„ e ^'"'^'■■■'"'"^ = 0, ai + bi = Ci + di (70) 



but also 



[ t; 7^ 7^ 1 dxi . . .dx„ f{S(xi, . . . ,Xn)) = 0, ai+bi=Ci + di (71) 

\dSai,...,a^ OSbi,...,b„ OSci,...,c„ OSdi,...,d„ J J V / 

with (almost) arbitrary choice of function f{S). Again, the validity of these differential equations obviously 
follows from the fact, that differential operator in the left hand side annihilates the integrand. As one can see, 
two different integrals - with integrands and f{S) - satisfy one and the same Ward identities. By itself, 
this fact does not guarantee equality of these two integrals. However, both integrals are SL{n) invariant and 
have one and the same degree of homogeneity: for any f{S) we have 



J dXl . . . dXn f(^XS(xi, . . ■,Xn)j = A "Z'' y dXl . . . dXn f(^S{xi, . . .,Xn)^ 



(72) 



Together, these three properties - common Ward identities, common degree of homogeneity and SL(n) invariance 
- can be considered as a proof (or, at least, as a strong evidence in favor) of equivalence of these two differently 
looking integrals: for arbitrary good function f{S), we have 



j f (^S{xi, . . . ,Xn)^ dxi . . . dx„ = const{f) J e ^^'''^'"'''^'^^ dxi . . .dx„ 



(73) 



The constant here depends on the choice of /, but does not depend on 5*. Property (73) is sometimes called 
"action-independence" of integral discriminants. Even before specifying the class of good functions f{S) and 
proving the action independence, let us consider a simple illustration with f{S) = e~^^: 

J J VAi ...Xn J VdetS 

/e-(«--'^^)' d"x= /e-(^-?+-+^-^)' , ^ f-U+...+.lY const' 

J J VAi ...Xn J VdetS 
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As one can see, these two differently-looking integrals are indeed essentially one and the same (they arc pro- 
portional with ^-independent constant of proportionality). To specify the class of good functions and to prove 
(73), let us make a change of integration variables 



Xl = p, X2 = pZ2, X3 = PZ3, ■ ■ ■ ,X„ = pZn 



i.e. pass from homogeneous coordinates Xi to non-homogeneous coordinates Zi = Xi/xi. Then 

J g-s(xi,x2,...,^„) dxi...dxn = J p"'^dp J dz2...dzr, e-"^^'^'"''^-'''^^ = 



„-l -p"- , \ f dZ2...dZn 

' ' '"'-y g(M.,. ..,.„)- 



/r 



For the right hand side of (73) we have 

j f(^S{xi,X2,---,Xn)j dxi...dx„ = J p^'^^dp J dZ2...dZn f(^P^S{l,Z2,...,Z„)j 



,Z2,...,ZnY''' 

Both integrals over p are just S-independent constants. As one can see, (73) is valid, iff these integrals 

j p^-^e-'^dp and j p''-^!{jf^dp 

are finite over one and the same contour. This condition specifies the class of good functions f{S). For example, 
all functions f{S) = exp ( — S'^) for > fall into this class. 

5.3 Evaluation of areas 

Action independence of integral discriminants allows to make different choices for f{S). This explains why 
integral discriminants appear in a variety of differently looking problems. For example, we can take 



m = e{i -s) = 



1, s<i 

0, S>1 



(74) 



which makes sence only in the real setting, i.e, when the contour of integration is real (and the form S is positive 
definite). This choice of f{S) provides an important geometric interpretation of integral discriminants: 



J g sixi,x2,...,x„) dg.^ dy.^ (60)^(61) j ^ volume, bounded by the hypersurface S = 1 



S<1 



(75) 
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This relation to bounded volume is a nice practical application of integral discriminant theory. As an example, 
let us consider the case of J214. Take a particular algebraic curve 



S{x,y) = {x +y) +txy =1 



(76) 



The area, bounded by this curve, by definition can be calculated in polar coordinates as 



\/l + e cos^ </> 



where K is the complete elliptic integral of the first kind [17] and can be represented as 



27r ^ / 1/2,1/2 



4 / e 



e + 4 



(77) 



From the integral discriminant side, the homogeneous quartic form S{x, y) has two invariants 



J 8 2 1 2 16 2 1 2 1 3 



(78) 



and the two branches of the integral discriminant are 



-1/4 



2F1 



1/12.5/12 
1/2 



(6 + 2)^(6 + 8)^(6-4)^ 

(6^ + 4e + 16)3 



(79) 



16 2 1 2 1 3 

— e e e 

9 3 6 36 



8 2 1 2 
3+r+6' 



-7/4 



7/12, 11/12 
3/2 



(e + 2)^(6 + 8)^(6 - 4)^ 
(e2 + 4e + 16)3 



(80) 



As noticed in [5] , only one linear combination of these two branches stays regular at e — > 0: namely. 



■^214^-^6 r(7/i2)r(ii/i2);"^2|4(e) lb 



+1/4 r(i/2) \ (2), , 
r(i/i2)r(5/i2) ) "^^^^ ' 



(81) 



Therefore, it is this branch that should be identified with the bounded area (regularity of the bounded area 
in the limit 6 — > is fairly obvious). Comparing the bounded area with the regular branch of the integral 
discriminant, we find a hypergeometric identity 



y(e) = const ■J^[f^(e) 



(82) 
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Eq. (82) ia indeed a valid hypergeometric identity, which can be seen by expanding both sides in powers of e. 



. 1^92 



25 a 1225 4 
7re + ..,„.„„. 7re + 



4194304 



(83) 



4[r'(^) 



i 9 2 

4 64*^ 4096^ 



25 

e + 



65536 



1225 4 

:£ + • 



16777216 



(84) 



what fixes const = 47r in (82). As one can see, the bounded area and corresponding integral discriminant indeed 
coincide, up to a constant factor. 

The above calculation with invariants is certainly not the most practical way to evaluate the integral dis- 
criminant: it is rather a conceptual demonstration of equivalence between these quantities, emphasizing the 
importance of correct choice of the regular branch. Usually, a much more convenient way to do is to perform a 
series expansion directly in the integral (66). Doing so automatically guarantees that a correct branch is chosen, 
and often allows to calculate the bounded volume in a very convenient way. As an illustration, consider the 
algebraic hypersurface 



a{xi + xi + xlY + b(xi + xl + x\)xi + cxi = H 



(85) 



in four-dimensional space. Using the relation with integral discriminants, we find the bounded volume 



hOO +00 +CXD +CXD 



a,h,c)= J dxi J dx2 J dxs J dxi exp ^-^{xl + xl + xl)^ - ^{xl + xl + xi)xi - ^xi 



-00 —00 —00 —00 



(86) 



Passing to new integration variables 



xi = pcos{6i) 
X2 = psm.{6\) cos(0) 
X3 = psin(^i) sin(0) 

XA=t 



with 



O<6i<7r,O<0<27r,O<p<oo, -00 < t < 00 



(87) 



and Jacobian sm{d), we obtain 
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oc 00 



V{a,b,c) ~ J J J p^dp J dtexp ^— — ^P^^^ ~ 



(88) 



The angular integration is easily done, and we find 



OC 00 



V{a,b,c) = 4:n j p'dp j dte^^ (^-'^p^ - ^p^t^ - ^t''] (89) 



-00 



Expanding into series 

^ OC 00 

y (a, b,c)=4nJ2^f p'+'' exp (- ^p^) dp f exp (- ^t') dt (90) 
'==0 ' -00 

and evaluating both integrals, we find 

, X &*= fk 1\ (k 3\ /a\-.V4-fc/2 /c\-l/4-fe/2 

n-'^'-) = 2 ^ + 4) ^2 + 4) (ft) in) (91) 

fe=0 



The sums are easily evaluated in elementary functions, and we finally obtain 

V{a, b,c) = ^ , ^ (92) 

In this example, the use of integral discriminant allows to obtain a simple formula for the volume in radicals. 
Thus integral discriminants are useful in application to the calculation of bounded volume. 

5.4 Algebraic numbers as hypergeometric series 

A characteristic property of integral discriminants, which is evident from definition, is SL(n) invariance; integral 



dx\ . . . dxn e 



-S(xi,...,Xn) 



does not depend on the choice of basis in the n-dimensional linear space and, consequently, depends only on 
SL{n) invariants Ik{S) (this dependence is made explicit in the above summary table). This property is not 
generic: many other quantities, which belong to the class of partition functions, are not SL{n) invariant. 
An important example is provided by roots of polynomials, also known as algebraic numbers: 



/o + /iA + ... + /.A'- = 



42 



This equation defines A as a function of fo,...,fr- 



A = A(/o,. ..,/.) 



and this function is naturally multi- valued - there are exactly r branches of A, since algebraic equation of degree 
r has exactly r roots in the complex plane. For low r this function can be expressed in radicals. For r = 2 this 
is merely the solution of quadratic equation: 



A(/0,/l,/2) 



1 /i ± V/f - 4/2/0 

2 /2 



Similarly, for r = 3 one of the roots is given by 



^/f f f 1 h , (36/1/2/3 - IO8/0/I - 8/1 + 12/3^12/^/3 - 3/f /I - 54/1/2/3/0 + 81/g/l + 12/o/|)^/^ 

3 /3 0/3 



6/1/3 - 2/1 



3/3(86/1/2/3 - IO8/0/I - 8/1 + 12/3^12/^/3 - 3/f /I - 54/1/2/3/0 + 81/2/1 + i2/o/|)i/3 

and the others are analogously expressed through CTibic roots. For r = 4, radicals of order 1/4 rtq already 
necessary, and for generic r expression in terms of radicals is impossible (as rigorously shown by Abel and 
Galois). Still, despite it is impossible to express them through elementary functions, for any r the algebraic 
numbers satisfy the same differential equations 





Ofa Ofb dfc Ofd ' 



(93) 



as integral discriminants do. For example, it is possible to check the Ward identities directly for r = 2: 



\ I fi ± Vf! - 4/2/0 



dfo df2 dfl dfl 



= 



similarly for r = 3 



'_d__d 

.dfo 0/2 Oh Oh 



1 h {36hhh - IO8/0/I - 8/1 + 12/3712/3/3 - 3/f /I - 54/1/2/3/0 + 81/oVI + 12/o/|)V^ 
'3/3 6/3 



6/1/3 - 2/1 



= 



3/3(36/1/2/3 - IO8/0/I - 8/1 + 12/3 Vl2/?/3 - 3/iVI - 54/1/2/3/0 + 81/^/1 + 12/o/|)i 



/3 
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d_Jl d_jr\ ( 1 /2 ^ (36/1/2/3 - IO8/0/I - 8/1 + 12/3^12/^/3 - 3/iVI - 54/1/2/3/0 + 81/g/l + 12/o/|)^/^ 



dfodh dfi 9/2 M 3/3 6/ 



3 



6/1/3 - 2/1 



3/3(35/1/2/3 - IO8/0/I - 8/1 + 12/3^12/1^/3 - 3fifi - 54/1/2/3/0 + 81/^/1 + 12/o/|)V3 

and also for r = 4. For arbitrary r (when no expressions through radicals are available) validity of Ward identities 
(93) becomes obvious, if one recalls the conventional Cauchy integral representation for A: 

\{fo,---,fr)= ^ Z d\og(^fo+flZ+...+ frZ"-^ = - ^ log (/o + /l^ + . . . + frZ^^^ dz (94) 



All this - multivaluedness, Ward identities and integral representations - make algebraic numbers quite similar 
to integral discriminants. Of course, since roots are not SL-invariant, they are not literally equal to inte- 
gral discriminants. Still, as a consequence of Ward identities, solutions of algebraic equations are given by 
hypergeometric functions, similarly to integral discriminants. The basic example here is the case of r = 2: 



/o+/iA + /2A^ =0 (95) 



This equation has two independent solutions, given by explicit formulas 



\ (f f f \ - 1 /i + V/i - 4/2/0 X ( f f f \ _ 1 h - Vfi - 4/2/0 

-^+(/o, /i,/2j 7 , -^-(/o, /i,/2j 7 (yoj 

2/2 2/2 



If viewed as a series, these functions take a form 



fo /0/2 2/g/l 5/oVI 14/oVI 42/oVI 

'~ /i f? n n n ^ ' 



and 



/i ^ /o ^ fSf2 ^ 2/^/1 5/oVI ^ I4/0V2" ^ 42/g/l 



Given these series, one can also verify the validity of Ward identities term-by-term. Note again the similarity 
with integral discriminants: there are two branches, which satisfy one and the same differential equation and 
differ only by asymptotics (at large /i, one branch grows asymptotically as — /1//2, while the other tends to 
zero). The coefficients here are famous Catalan numbers, 
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1,2,5, 14,42,...= 



(2k)\ 



A;!(A; + 1)! 



(99) 



They axe integers and have a celebrated combinatorial meaning of counting the number of tree diagrams, which 
appear in the Bogolubov-like iteration procedure [16]. However, we are interested not so much in diagrammatical 
and/or combinatorial meaning of these numbers, but rather in their hypergeometric nature: the solutions are 
given by explicit hypergeometric series 



{2k) 



/2 



and 



/o p 

Jl 



4/2/0 



(100) 



J ^ I — n 



(2fc)! /o^+VI 

2fe+l 



/i , /o p 

J2 Jl 



J2] 



4/2/0 



(101) 



which just occasionally reduce to square roots. Note again the similarity with integral discriminants: singularity 
of these series is situated at the unit value of argument, when /i = 4/2/0, i.e, when the algebraic discriminant 
vanishes. In variance with expressions through radicals, expressions through hypergeometric functions can be 
written for arbitrary r: say, an equation + x + c = which is not solvable in radicals, can be solved in series 
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Looking at these numbers, it is easy to find 



1, 6, 51, 506, 5481, 62832, 749398, . . . = 



(6fc)! 



fc!(5fc + 1)! 



(103) 



These numbers are again integer and enumerate sextic (6-ary) trees (rooted, ordered, incomplete) with k vertices 
(including the root). Moreover, their generating function is again hypergeometric: 
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Note, that this function represents only one of the roots: the other roots are represented by other hypergeo- 
metric functions, which are completely analogous and which we do not include here. For more details about 
hypergeometricity of solutions of algebraic equations, see [18]. Note also, that integerness of the coefficients is 
not a generic property, it is due to particular form of the above examples. Say, algebraic equation 



(a^ + o)^ + (a; + 6)* = 
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has a solution 



x = 1+ ^-a + 2b- %^o^ - 606 + 376^ + . 
2 8 



with already non-integer coefficients and exact formula 
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Let us make one more short note about SL{n) invariance. As aheady mentioned, the roots theirselves are 
not SL{2) invariant: under transformation G € SL{2) they transform as 

^ ^ G11G22 - G12G21 = 1 (105) 

1^21 A + Lr22 



Thus they are not expressible through the elementary invariants. However, the double ratios 



(\i - Aj)(Afc - A;) 

= (A.-A.)(A,-AO 



are invariant under SL{n) transformations (this is easy to check directly). For example, for r = 4 

/o + /iA + /2A' + /3A^ + /4A* = (107) 

there are two elementary invariants 

h = 2/0/4 - i/1/3 + i/l (108) 



h = /0/2/4 - pofi - lf!f4 + i/1/2/3 - ^/l (109) 



and a single double ratio A1234, which is SL(2) invariant and thus expressible through I2 and 73. To express 
A1234 through the elementary invariants, let us consider a particular case 



/o = -A", /i = e, h = 0, /3 = 0, /4 = 1 (110) 
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For these values of parameters h 

Ai = 



-2A^ and h = — 3e^/8. It is easy to obtain, that 
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These four branches correspond to e-deformations of the four roots of the simple equation A* = A*. Accordingly, 
the double ratio can be calculated as a series in e and equals 



_ (Ai - A2)(A3 - A4) _ 1 3/ 337 
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Expressed through invariants, it takes form 
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This is the simplest way to express the double ratio through h and h. Because of its SL{2) invariance, this 
equality will remain the same for arbitrary /o, . . . , /4, not just for particular values (110). 

Let us finish this section by noting, that all algebraic numbers - in other words, the set of numbers which are 
solutions of some algebraic equations - form a field. This in particular means that the set of all algebraic numbers 
is closed under addition and multiplication. For example, suppose that xi and X2 are algebraic numbers. I.e, 
suppose that xi is a root of f{z) and X2 is a root of g{z). Then polynomials f{z) and g{t — z) have a common 
root for t ^ xi+ X2, so that xi + 0:2 is a root of resultant of f{z) and g{t — z). Since resultant is a polynomial, 
xi + X2 is an algebraic number. Similarly one can prove that xi — X2, x\X2 and xi/x2 for X2 ^ are algebraic 
numbers. Thus, algebraic numbers actually form a field. 

By analogy, one could in principle consider 2-algebraic numbers {x,y): solutions of systems of equations 




where / and g are some non-homogeneous polynomials, and their straightforward generalisations for arbitrary 
number n of independent variables. However, multiplication and division of such multi-component objects 
seems to be unnatural. A more natural direction of generalisation makes use of the similarity between algebraic 
numbers and integral discriminants: it is interesting to understand, in what exact sence the set of integral 
discriminants of homogeneous forms is a ring or even a field. This brings us directly into the field of motivic 
theories (see [19] for a recent review "for physicists") which is, however, beyond the scope of this paper. 
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6 Conclusion 



We have briefly described (hopefully in explicit and elementary way) several research lines in resultant theory. 
These include Sylvester-Bezout (determinantal) and Schur (analytic) formulas for the resultant, invariant formu- 
las and symmetry reductions for the discriminant, Ward identities and hypergeometric series representations for 
partition functions (integral discriminants and roots of algebraic equations). Unfortunately, several important 
directions - such as Macaulay matrices [10], Newton polytopes [3], Mandelbrot sets [20] and others - remain 
untouched in this review. They will be discussed elsewhere. 

Let us emphasise that the objects and relations, which we consider in present paper, are somewhat under- 
investigated. Many interesting properties of resultants, discriminants and especially of partition functions have 
attracted attention only recently. Moreover, new unexpected properties and relations continue to appear. All 
this indicates, that resultant theory and the entire non-linear algebra is only in its beginning. 
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